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Abstract 

Some integral properties of Jack polynomials, hypergeometric functions and invari- 
ant polynomials are studied for real normed division algebras. 



1 Introduction 
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During the 1 960s, real and coni plex zonal polynom ials were studied exhaustively bv I Jame j 
()l96ll . ll964 , IConstantind ('1963') and'Khatri' ('1966'), among man y others . Excellent reference 
books include those bv [Muirhcad (1982), Takemura_ (.198j), iFarrelll ([1985, ') and Mathail 
(|l997l ). which summarise many of the results published to date. 

Hypergeometric functions with a mat rix argument were fi rst studied bv lHer j (1955) and 
defined in terms of zonal polynomials bv lConstantine ( 19631 ). Hypergeometric functions of 
one or two matrix arguments have been applied in many are as of science a nd technology, 
including multivariate statistical analysis (^Muir head (1982) and lMatha j (199^)), ran dom ma- 
trix theory (Mctha (1991) and Forrester (2009) ), wir eless communication s (Ratnar aiah et all 
(l2005a|y) shap e theo r y iGoodall and Mardial (1 1991 and ICaro-Loper~et al\ (l2009h) . 

Later iDavid (|l979l ). iDavid (|l980f) . IChikusd (|l980l ) and lChikuse and Davislll986l ). intro- 
duced a class of homogeneous invariant polynomials with two or more matrix arguments, 
which generalise the zonal polynomials; many of their basic and integral properties are 
studied in real cases. 

In the context of multivariate statistics, zonal polynomials were initially used to express 
many noncentral matrix variate distributions. However, there were other distributional 
problems that could not be solved using zonal polynomials. In these latter cases, invariant 
polynomials were used to obtain explicit expressions of doubly noncentral matrix variate 
distributions, matrix variate distribution functions an d the jo i nt de nsity of eigenvalue s of 
matrix variate beta type I and II distributions, etc. see IJamesI (|l964[) and lDavisI (|l980l) . 

During the 1990s, zonal polynomials regained prominence but from a more general point 
of view, in which it was observed that zonal po l ynom ial fo r real and complex cases ar e 
particular cases of Jack polynomials, see ISawyeil (1992) and Goulden and JacksonI ( 1996f l. 
among many others. In terms of Jack polynomials, it is possible to give a general definition 
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for hypergeometric functions, see I Gross and Richards! (|l987l) . iGross and Richards! (|l989! ). 



among many others. In both cases, Jack polynomials and generalised hypergeometric func- 
tions are written in term of a parameter, denoted by a or /3 and with which, for example 
for /3 = 1, 2 o 4, the zonal polynomials and hyp ergeometric function s are o b tained for real 
compl ex a nd quaternion cases, respect ively, see ICross and Richards! (1987), Dimitriu et al 
(|2005l ) and lKoev and Edelman ! (|2006!) . 

The properties for Jack polynomials and hypergeometric functions with a matrix argu - 
men t have b een studied by iHerj (|l955t) . jjamel (|l964!) . [Constantind (Il963!). iKhatril (Il966!) 
and iMuirhea d (1982) 
(Il984l ). iFarrelll (Il98 



zona l polyno mials); by iJames (Il964ir iTakemura 



in t he real case . ^ ^ ^ ^ 

^ , , — and Ratnaraiah et al. ( 2005al lbl) in the complex case (Schur func- 
tions); by Li and Xud (2009 ) in the quaternion case and bv I Gross and Richardd (1987) and 
Caro-Lopera et aZT ( 2007t ) in the general case (real, complex and quaternion cases), among 



many others. 

A serious obstacle encountered when Jack polynomials, hypergeometric functions and 
invariant polynomials are to be use d is the question of their ca lcula tion. Fortun ately, with 
the excellent algorithm proposed by Koev and Edelman ( 2006!) and Koev (l2004l). it is now 
possible to use these tec hniques in many applications, see Ratnaraiah et a/] ~( 2005b! ) and 



Caro-Lopera et al. 



()2009l) . Unfortunately, this obstacle remains for the general case of in- 
variant polynomials. 

Let us take into account that there are exactly four normed division algebras: the real 
numbers (K), complex numbers (£), quaternions (Sj), and octonions (D); moreover, these 
are the only al ternat i ve div ision algebras, and all division algebras have a real dimension o f 
1,2,4 or 8, see lBaezI (|2002L Theorems 1, 2 and 3). Furthermore, according to iBae d (|2002l ). 
there is still no proof that octonions are useful for understanding the real world. 

In this paper, we generalise diverse integral properties of Jack polynomials, hypergeomet- 
ric functions and invariant polynomials for normed division algebras. Note that we can only 
conjecture the results for the octonion case, because many o f its related matrix problems 
are still under study. However, for example in [Forrester ()2009L Section 1.4.5, pp. 22-24) it is 
proved that the bi-dimensional density function of the eigenvalue, for a 2 x 2 octonionic ma- 
trix with symmetric normal distribution, is obtained from the general joint density function 
of the eigenvalues for the symmetric normal distribution, assuming to = 2 and /3 = 8. The 
material in the present paper is organised as follows: Section [2] provides some notation and 
preliminary results about Jacobians, gamma and beta multivariate functions and invariant 
measures. The definition and many integral properties of Jack polynomials are obtained in 
Section [3l Many extensions of the integral properties of hypergeometric functions with one 
and two arguments are studied in Section [H For invariant polynomials with two matrix 
arguments, in Section [5] we derive diverse integral properties, such as the inverse Laplace 
transformation, gamma and beta integrals, etc. Finally, in Section [SI we show diverse ap- 
plications of some results derived previously, such as the distribution function of a central 
Wishart distribution for normed division algebras, its joint eigenvalue density and the dis- 
tribution function of the largest and smallest eigenvalues. We emphasise the conditions that 
must be met by the parameters that take part in many integral properties in the cases dis- 
cussed, because, even in the original references, these conditions were omitted or established 
inexactly. 



2 Preliminary 

Let Cf^ „ be the linear space of all n x to matrices of rank m < n over ^ with to distinct 
positive singular values, where ^ denotes a real finite- dimensional normed division algebra. 
In particular, let GL{m,^) be the space of all invertible m x m matrices over ^. Let 
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be the set of all n x m matrices over ^. The dimension of over 3? is /3mn. And let 

us recall that the parameter /3 has traditionally been used to count the real dimension of 
the underlying nor med division algebra. In other branches of mathematics, the parameter 
a = 2/(3 is used, see lBaezj ()2002l ). 



Table 1: Values of the /? = 2/a parameter. 





a 


Normed divison algebra 


1 


2 


real (5R) 


2 


1 


complex {€.) 


4 


1/2 


quaternion {?)) 


8 


1/4 


octonion (D) 



T 

Let A G 5'"^™, then A* — A denotes the usual conjugate transpose. The set of matrices 
Hi G such that H*Hi = is a manifold denoted „, termed the Stiefel manifold 

(Hi are also known as semi- orthogonal {(3 = 1), semi- unitary (0 = 2), semi-sym plectic 
(/3 = 4) and semi- exceptional type {f3 = 8) matrices, see iDrav and Manogud (|l999l )V The 



dimension of „ over K is [f]mn — m{m — l)/3/2 — to]. In particular, „ with dimension 
over 5ft, [m{m+l)(3/2~m], is the maximal compact subgroup il'^ (to) of C!^ „j and consists of 
all matrices H e grmx™ such that H*H = 1^. Therefore, ii^{m) is the real orthogonal group 
0{m) {P = 1), the unitary group lA{m) {j3 ~ 2), the compact symplectic group Sp{m) (/3 = 4) 
or exceptional type matrices Oo{m) (/3 = 8), for 5^ = 5ft, £, ^ or D, respectively. Denote by 
the real vector space of all S G ^"^^"'^ such that S = S*. Let be the cone of positive 
definite matrices S G g^^^™; then is an open subset of 6^. Over 5R, 6^ consist of 
symmetric matrices; over €, Hermitian matrices; over ij, quaternionic Hermitian matrices 
(also termed self-dual matrices) and over D, octonionic Hermitian matrices. Generically, 
the elements of are termed as Hermitian matrices, irrespective of the nature of g". 
Thcd imcnsion of over 5ft is [to (to — 1)/3^"2]/2. Let be the diagonal subgroup of 

consisting of all D G D — diag(c?i, . . . ,dm)- Let 'Z^ijn) be the subgroup of all upper 

triangular matrices T G such that ty = for 1 < z < j < to; and let 1^{m) be the 

opposed lower triangular subgroup T^(to) — ^T^(to)^ . For any matrix X G :J"x™j rfX 

denotes the matrix of differentials (dxij). Fi nally, we define t he measure or volume element 
(dX) when X G S-™'^",©^, S)^ or V^„, see iDimitriul (l2002l) . 



If X G g'"^™ then (dX) (the Lebesgue measure in g'"^*™) denotes the exterior product 
of the f3mn functionally independent variables 



(fe) 

1 j=l fe=l 



(dx) = A A A 

Remark 2.1. Note that for G 5^ 

da;^- = A '^^ 
In particular for = 5ft, £, f) or D we have 





(fe) 
y ■ 

k=l 
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• Xij G then 



dxi 



/\ dx 



(k) 



(2) 

ix)/ e £. then 



(fe) 



J2) 



• (3) 



dxi 



- A:Xj-^' e io, then 



: dx[f A dxg^ A dx,g^ A dx[f = /\ dx|f. 



fc=i 



J2) 



„(3) 



• Xy = x^y + eix\y + e2xY/ + eaxl^'' + e4X-,-' + esx- + eex 



(5) 



(6) 



(7) 



■ ejxlj £ D, then 



dxi 



dx,^,^^ A dxlf A dx,^f A dx|t^ A dx^f A dx^f A dx,^^ A dx^f ^ A dx 

''J ''J ''J ''J ''J ''J ''J ''J / \ 

k=l 



If S e 6^ (or S e then (dS) (the Lebes gue measure in 6^ or in 1.^{m)) denotes 

the exterior product of the m{m+ I) (3/2 functionally independent variables (or denotes the 
exterior product of the m(m — l)(3/2 + n functionally independent variables, if sa G for 
alH = 1, • ■ • , m) 

m /3 

(dS) = I '-"'=^ „ 

f\dsu f\ f\ds[';\ ifs,, g3?. 

i—l i<j k—1 

Remark 2.2. Since generally the context establishes the conditions on the elements of S, 
that is, if Sij G 5R, G £, G i3 or G D. It shall be considered 



m 13 



m P 



idS) = /\ f\ds[f^ /\dsu /\ f\ds 



(k) 



i<j k=l 



i<j k—1 



Observe, too, that for the Lebesgue measure (dS) defined thus, it is required that S G 
that is, S must be a non singular Hermitian matrix (Hermitian definite positive matrix). In 
the re a l case , when S is a positive semidefi n ite ma t rix, its corresponding measure is studied in 
Uhlid (Il994l). iDiaz-Garcia and Gutierrej (Il997l) . iDfaz-Garcfa and Gonzalez-FarTai (|2005ah 



and iDiaz-Garcia and Gonzalez- Fariad (|2005bl ) under different coordinate systems. 

If A G then (dA) (the Legesgue measure in 2)^) denotes the exterior product of the 
Pm functionally independent variables 



n 



(dA) = f\f\d\ 



(fe) 



i=l fe=l 



If Hi G „ then 



(HJdHi) 



A A 

i=i j=i+i 
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where H = (H1IH2) = (hi, . . . , h,„|h„i+i, . . . , h„) £ iX^{m). It can be proved that this 
differential form does not depend on the choice of the matrix H2 and that it is invariant 
under the transformations 

Hi^QHPi, Q andP eil'3(m). (2.1) 

When m = 1; vf „ defines the unit sphere in 5^". This is, of course, an {n — dimensional 
surface in S"". When m — n and denoting Hi by H, (H*(iH) is termed the Haar measure 
on il^ (m) and defines an invariant differential form of a unique measure i' on il^ (m) given 

by 

zy(M) = / (H*dH). 
Jm 

It is unique in the sense that any other invariant measure on il^ (m) is a finite multiple of 1/ 
and invariant because is invariant under left and right translations, that is 

The surface area or volume of the Stiefel manifold V,fj „ is 

Vol(VlJ = / (HldHi) = (2.2) 

and therefore 

(dHi) = 2^(HtdHi) = l|M^(HtdHi). 

is the normalised invariant measure on V,'^,„ and (dH), i.e. with ( TO = n), it defines 
the normalised Haar measure on il'^(TO). In (|2.2p . r^[a] denotes the multivariate Gamma 
function for the space 6^j, and is defined by 

Ti[a] = f etr{-A}|A|°-("-i)''/2-i(dA) 

m 

^ 7r"("-i)'5/4[|r[a- (i- l)/?/2] 
1=1 

m 

^ ^™(™-i)/3/4]^r[a- (to-i)/3/2], (2.3) 
1=1 

where etr|-| = exp|tr(-)j , | • | denotes the determinant and Re(a) > (to — l)/3/2, see 



wnere etr-j-j- = exp - itr(-)K | • | aenotes tne aetermmant ana tie[a) > [m — i)p/z, see 
Gross and Richards! (1983). This can be obtained as a particular case of the generalised 
gamma function of weight k for the space 6^^ with k = (fci, fc2, . . . , km), ki > > ■ ■ ■ > 
km > O1 taking n = (0, 0, . . , 0) and which for Re(a) > (to — l)/3/2 — km is defined by, see 
Gross and Richards! |l987l ). 



r^[a,«] = / etr{-A}|Ar-('"-i)'3/2-i^^(A)(dA) 

m 

^ TT^^^-i^'^/^J^ria + fc, -(i-l)/3/2] 



i=l 



(2.4) 
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where for A 6 &fL 



(2.5) 



i=i 



with Ap = (ors): f^s = 1, 2, . . . ,p, p = 1, 2, . . . , m is termed the highest weight vector, see 
Gross and Richardd(|l987l ). 



Remark 2.3. Let V{&^) denote the algebra of all polynomial functions on 6^, and Vki&m) 
the subspace of homogeneous polynomials of degree k and let V^{&!^) be an irreducible 
subspace of 7^(6 such that 

Pfe(6f„) = ^07'«(6^). 

Note t hat q„ is a homogeneous polynomial of degree fc, moreover q„ e 'P''((3^), see Gross and Richardi 
(Il987h . 

In (|2.4p . [a]^ denotes the generalised Pochhammer symbol of weight k, defined as 



^m(m-l)/3/4 "Q p[„ _^ _ ( • _ ;^)^/2] 
1=1 



where Re(a) > (m — l)/3/2 — fc,„ and 

{a)i = a{a + 1) • • • (a + i - 1), 

is the standard Pochhammer symbol. 

A variant of the generalised gamma function of weight k is obtained from iKhatril ( 19661 ) 
and is defined as 



r/3 



[a,-K] = [ ctr{-A}|Ar-('"-i)'3/2-ig,(A-i)(dA) 



nr[a-A:, - (m-z)/3/2] 



i=l 



^ 7r™('"-i)/3/4j^r[a-/c, -(^-l)/3/2] 

i=l 



'a+{m-l)f3/2 + l]Z' 



(2.6) 



where Re(a) > (m - l)/3/2 + fci. 

The two expressions of r^[a, ], r(^[a, k] and r^[a, — k] as the product of ordinary gamma 
functions are obtained using the proofs corresponding to A = TT* and A = T*T with the 
corresponding Jacobian given in Lemma l2.5l Alternatively, note that for any function g{y) 



Y[g{x + i - 1) = Yldix + q-i), 



(2.7) 



and 



Y[g{x - I + 1) =Y[g{x - q + i). 



(2., 
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Similarly, from Herd ( 1955 . p. 480) the multivariate beta function for the space 6^^, can 
be defined as 



0<S<I,„ 

[ |R|"-(™-i)''/2-i|i^ ^ R|-('^+'''(dR) 
ri[a]Ti[b] 

ri[a + b] ' 



(2.9) 



where R = (I - S)-i - I, Re(a) > (m - l)/3/2 and Re(6) > (m - l)f3/2 

Some Jacobians in the quaternionic case are obtained in iLi and Xuel(l2009l). We now cite 
some Jacobians in terms of the parameter /3, based on the work o fiDimitriul diooi). We also 
include a parameter count (or number of functionally independent variables, #fiv), that is, 
if A is factorised as A = BC, t hen th e parameter count is written as #fiv in A = [#fiv in 
B] + [#fiv in C], see lDimitriul (|2002[ ). 



Lemma 2.4. Let X and Y £ 'P^, and let Y = AXA* + C, where A and C G -C^.m are 
constant matrices. Then 

(dY) ^ |A*A|'5("-i)/2+i(rfx). (2.10) 

Lemma 2.5 (Cholesky's decomposition). Let S G and T G 1,^(m) with tu > 0, 
i = 1, 2, . . . , TO. Then 

• parameter count: (3m{m — l)/2 + ?7i = [3m{m — l)/2 + m and 



(dS) 



2"'Y[t 

1=1 

m 

2" 



/3(m-i) + l 



(dT) if S = T*T; 



/3(i-i) + l 



(2.11) 



(dT) if TT* 



Lemma 2.6 ( Spectral decomposition). Let S G ^^j- Then, the spectral decomposition 
can be written as S = WAW*, where W G il''(m) and A — diag(Ai, . . . , A™) G S)™, lyzt/i 
Ai > • • • > Am > 0. Then 

• parameter count: Pm{m — l)/2 + to = [/3TO(m + l)/2 — to — (/3 — 1)to] + [to] and 

m 

{dS) = 2-"'t^b J^(^x, - Xjf{dA){W*dW), (2.12) 



where 



0, /3=1 

-TO, /3 = 2 

-2to, /3 = 4 

-4to, l3 = 8. 



3 Integral properties of Jack polynomials 

In this section we review and study diverse integral properties of Jack polynomials for normed 
division algebras. However, let us first consider the following remarks and definitions. 
Remark 3.1. Note that Jack polynomial s and hyp ergeometric functions with one or two 
matrix arguments are valid for /3 > (jKoev and E dclman (2006 )), but in our case (3 
denotes the real dimension of ^. Also, we use the parameter /3 indeed of a in the definition 
of the Jack polynomials and hypergeometric functions, with the equivalence shown in Table 

m 
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Then: 

Let us cha racterise the Jack symmetric function jjf\Xi, . . . , Am) of parameter f3, see 



I Saw veil ( 1997 ). A decreasing sequence of nonnegative integers k = (fci,A:2, . . .) with only 



finitely many nonzero terms is said to be a partition of fc = ^ /Cj. Let k and t — (ti,t2, ■ ■ ■) 
be two partitions of k. We write r < k if J2l=i — J2i=i ^'^^ each t. The conjugate of k is 
k' = {k[,k2, ■ ■ ■) where k'^ = card{j : kj > i}. The length of k is l{k) = max{i : fc^ 7^ 0} = k[. 
If 1{k) < m, it is often written that k = (/ci, fe, • • ■ , fcm)- 

The monomial symmetric function Mk(-) indexed by a partition k can be regarded 
as a function of an arbitrary number of variables such that all but a finite number are 
equal to 0: if Ai = for i > m > 1{k) then Mk(Ai, . . . , A^) = ^ A^^ • • • Af^, where 
the sum is over all distinct permutations {61, . . . ,S„i} of {ki, . . . ,km}, and if 1{k) > m 
then Mk(Ai, . . . , A„i) = 0. A symmetric function / is a linear combination of monomial 
symmetric functions. If / is a symmetric function then /(Ai, . . . , Am, 0) — /(Ai, . . . , Am)- 
For each to > 1, /(Ai, . . . , Am) is a symmetric polynomial in to variables. 

Then the Jack symmetric function jjf\Xi,...,Xm) with a parameter /3, satisfies the 
following conditions: 

jf)(Ai,...,Am) = ^i.„,.(/3)Af,(Ai,...,Am), (3.1) 
4^\^,---,l) = [-^j l[{{m~^ + l)^3/2),^, (3.2) 
2?^j('5)(Ai,...,Am) = ^fc,(fc.-l + /?(m-i))j(«(Ai,...,Am). (3.3) 

i=l 



where 



™ 52 ™ Id 

2 = E + ^ aT^ 9A" 



i=l i=l 



Here, the constants Vk,.t{I3) do not dependent on X[s but on k and r. Note that if 
TO < 1{k) then J,l'^''(Ai, . . . , Am) = 0. The conditions include the case /? = and then 
J,l°-'(Ai, . . . , Am) = e^' n" i("^ ~ * + l)*"': where e«;(Ai, . . . , Am) = n'=i efe,(Ai, . . . , Am) are 
the elementary symmetric functions indexed by partitions k, if to > 1( k) then e r (Ai, . . . , Am) = 
Sii<j2< - <v Ail • • • 1 and if TO < 1{k) t hen er(Ai, . . . , Am) = 0, see ISawveil (|l997l ). 

Now, from lKoev and Edelman ( 2006f ). the Jack functions 



jf)(X) = j(«(Ai,...,Am), 

wher Ai, . . . , Am are the eigenvalues of the matrix X G &m^ can be normalised in such a way 
that 

^Cf(X) = (tr(X))^ (3.4) 

or equivalently, such that 

f:E^=etr{X}, (3.5) 

k=l K 

where (X) denotes the Jack polynomials (for simplicity, we have replaced {0) by (3 as the 
superindex for the Jack polynomials) . These are related to the Jack functions by 



Cf(X) = |J|j(«(X), (3.6) 
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where 



(i,j)eK, 



and h1{i,i) = kj — i + 2{ki — j + 1) / (3 and h'^{i,j) — kj — i + 1 + 2{ki — ])/ (3 are the upper 
and lower hook lengths at £ k, respectively. Also, observe that for X = S*S and 

Y = W*W we have 

Cf (WXW*) = (SYS*). (3.7) 

, 2 



In particular for A^/^ such that (A^/^)~ = A 

(Y1/2XY1/2) = C^(Xi/2YXi/2). 



(3.8) 



Therefore, given that XY, YX, Y^/^XY^/^ and X^/^YX^/^ all have the same eigenvalues, 
we opt for convenience of notation rather than strict adherence to rigor, and write (XY) 
or Cf (YX) rather than Cf (Y^/^XY^^^), even though XY or YX need not lie in S^. Note 
that 

(Zi/^XZ^/^Y) = Cf (XZ1/2yz1/2), (3.9) 

for all X,Y e S^, and Z e From i Gross and Richard^ (|l987l Equation 4.8(2) and 

Definition 5.3) we have 



Cf(X) = Cf(I) 



HGU'3(m) 



(zJH*XH)(dH) 



(3.10) 



for all X G where (dH) is the normalised Haar measure on i}fi{m). Finally, for the c 
constant we have that C;^(cX) = c'^Cf (X). 

Some basic integral properties are cited below. For this purpose, we utilise the complex- 
ification &Pf = + i6f„ of 6f„. That is, consist of all matrices X G (5-«:)™x™ of 

the form Z = X + iY, with X, Y G S^^. We refer to X = Re(Z) and Y = Im(Z) as the 
real and imaginary parts of Z, respectively. The generalised r ight half-plane $ = + 
in e^*^ consists of all Z G 6^*^ such that Re(Z) G see (Cross and Richards! |l987l p. 
801). 

For any X,Y G 6^, 



/ (XH*YH)(dH) = 



'HGil'3(m) 

For all R G 6f„, Z G * and Re(a) > (m - l)/3/2 - k 



(3.11) 



/ etr{-XZ}|X|"-("-i)'3/2-iC^(XR)(dX) 



r^^[a,^]|Z|-"Cf(RZ-i) 
[a]fr^[a]|Z|-'^Cf(RZ-i). 



(3.12) 



Remark 3.2. I n general, the result (13.121) has been e s tablished under the condi tion, Re(a) > 
(m-l )/3/2, see lConstanth^(|l963h . rMuirheadl(fl983 ). lRatnaraiah a/.l(|2005bh and fLi and Xud 
(l2009l ). but in reality the correct condition is Re(a) > (m — l)/3/2 — km- This fact is im- 
mediate, observing that [a]^r^[a] — r^[a,K] and the different expressions for F^Ja, k] in 
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Let Re(a) > (m - l)/3/2 and Re(6) > (m - l)/3/2. Then 

/■ |-j^|a-(m-l)/3/2-l|j _ X|''-(™- l^/'/^-l Cf (XR) (dX) 

Jo<x<i 

rl[a,K]rf„[6] 



[a + 5, k] 



[a + b] 







(3.13) 



for all R e ©m'^; see Gross and Richardi ( 1987 . Theorems 5.5 and 5.9 and Corollary 5.10) 
for real, complex and quaternion cases. 

Remark 3.3. Observe that result (j3.13l) was establishe d under the condit i ons Re(a) > (m — 
1)0/2 a nd Re(b) > (to - })0/2 see iConstantind (|l963f ). iMuirheadI (|l982h . lRatnaraiah et all 
()2005bl ) and lLi and x"u3 ()2009h . but the correct conditions are in fact Re(a) > (to — l)/3/2 — 
km and Re(6) > (to — l)/3/2. This fact is verified by observing that [a]fr(^[a] = r^[a,K] 
and the different expressions for Tf^[a, k] in p.4p . 

We now extend several integral properties of zonal polynomials in the real and complex 
cases to normed division algebras. O ur first result is a genera lisation of one studied by 
Teng et al. I (|l989l) for real case, see also ICaro-Lopera et all (|2009). From this result, we can 



obtain diverse particular integral properties of Jack polynomials. 

Theorem 3.4. Let Z e ^ and U e 6f„. Assume 7 = f{z)z''"'-''-'^dz < 00. Then 
f /(trXZ)|X|''-("-i)'3/2-i(^^ (X-^U) (dX) 



T[am — k] 

for Re(a) > (to - l)/3/2 + ki, and 

[ /(trXZ)|X|'^-(™"i)''/2-iC;f (XU) (dX) 



izr"c,^(uz).. 



(3.14) 



w/iere i? = /(z)z"™+''"Mz < 00, Re(a) > (to - l)/3/2 - fc™ and k = (fci 
= fci + • • ■ + fcm. 



(3.15) 

,km) and 



Proof. Denote the left side of (|3.14p by /(U, Z). By (|3.10p and interchange of order on 
integration 

J(I,I) = / /(trX)|X|"-('"-i)''/2-icf (X-i) (dX) 
= Cf(I) / /(trX)|Xr-(™-i)''/2-i 

xl f <z4H*X-iH)(dH) ) (dX) 

= Cf(I) / /(trX)|Xr-(™-i)'^/2-ig4X-i)(dX). 
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Let X = TT*, from Lemma 

m 

Then 

J(I,I)= 2^0^(1) f \^^tlj\{{tu?'-''^ 



.fc,-(m-i)/3/2)-l(-^rp^j^ 



Applying the Fang and Zhand ( IQQOl . Lemma 2.4.3, p. 51) we obtain 



1(1,1) = Cf(I)- 



^m(m-l)/3/4 "Q ^[a - k, - {m - i)P/2] 



T[am — k] 



■ 7 



= '-■KWTTr rT-7, 

i [am — k\ 

with 7 = f{z)z°-"^~^~^dz < oo. 

Next, since the function X(U, I) is invariant under il'^(m) and in V^{&f^), there exists 
a constant d such that X(U, I) = dC/f(U), U e 6^,. It is obvious that d = 1(1, 1)/Cf (I), 
then 

J(U,I) = ff-[°''"] cf(U).7. 
i [am — fcj 

Now, let Z G and make the change of variable X Z^^/^XZ^^^^ in the integral 
defining X(U, Z). Then by jSH) 

i(u,z) = |zr°i(zi/2uzi/2,i), 

and hence 

Z(U, Z) = Ipt||Z|-Cf (Z1/2UZ1/2) . ^. 

i [am — k\ 

Therefore, for Z e and U e 6^ 

I(U,Z) = l|>l^|Z|-Cf(UZ) .7. 
i [am — fcj 

The result in (I3.14p now follows by analytic continuation in Z from to # = ^P^j + «6(^. 
The result in (I3.15P is obtained in a similar way. □ 

Now, by definition if k = then [a]^ = 1 and C^(X) = 1 from where: 
Corollary 3.5. Let Z e $ G 6^^. Assume 7 = f{z)z°-"^^^dz < 00. Then 



f /(trXZ)|X|°-('"-i)^/2"i(dX) = 1|J4|Z|-'^ -7, 
7xe<Bg. r[amj 



(3.16) 



/or Re(a) > (to - l)(3/2. 

If we take f{y) = exp{— y} in Corollary [33] we obtain 



/ etr{-XZ}|X|'^-("-i)''/2-i(dX) =^^[a]|Z^^ (3.17) 

and if Z = I we obtain the multivariate gamma function for the space 6^. 
Other particular results of Theorem [33] are summarised below. 
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Corollary 3.6. Let Z e ^ and U G e^. 

/ etr{-XZ}|Xp-("'-i)^/2-icf (X^^U) (dX) 

= ri[a,-^]\Z\~-C^JVZ), (3.18) 

/orRe(a) > (m - l)/3/2 + fci . 

Proof. This is obtained by taking f{y) = exp{— y} in (|3.14p . □ 
Corollary 3.7. Let Z G $ and U G and j G suc/i that Re(ma + j — A:) > 0, i/ien 

/" etr{-XZ}(trXZp|X|"-("-i)^/2-iCf (X-^U) (dX) 

= ^"[°'-"]^["'" + ^''^] |Z|-Cf(UZ), (3.19) 

/or Re(a) > (m — l)/3/2 + ki. And if j is such that Re(ma + j + fc) > 0, then 

[ etr{-XZ}(trXZ)^ |Xp-("'^i)^/2-iCf (XU) (dX) 

= ^"["'i^["^"+/' + ^]|Z|-Cf(UZ-^), (3.20) 
r[ma + k\ 

for Re(a) > (to - l)/3/2 - /c^. 

Proof. The desired resuh is obtained by taking /(y) = exp{— in Theorem 13.41 □ 
Corollary 3.8. Let Z G * and U G 6^ and 77 > t/ien 

/ (l + 2ry-4rXZ)-'^(™+")|X|"-("-i)'^/2-icf (X-^U) (dX) 

|Z|-''Cf(UZ), (3.21) 



_ r^[a, -«:]r[(/3 - l)aTO + I3tj + k] 



(277-i)"™-'=r[/3(TOa + 77)] 
for Re(a) > (to - l)/3/2 + ki. And 

f (1 + 277-1 tr XZ)-'^(™+'') |x|'^-(™-l)/3/2-l(;;^ (XU) (dX) 

Jxeqjf^. 

Zr'^Cf(UZ-i), (3.22) 



'xeqjj^ 

r^[a,K]r[(/3-l)aTO + /377-fc] 



- (27;-i)'"™+'=r[/3(ma + ?7)] 
/or Re(a) > (to - l)/3/2 - /c„. 

Proof. The desired result is obtained by taking f{y) = (1 + 2r]^^y)~^'^'^"^^^^ in Theorem 

inn □ 



Many other interesting particular cases of Theorem 13.41 can be found, for example by 
defining /(trXZ) as t he k ernel of matrix variate generalised Wishart distributions, see 
Fang and Zhand (|l990[ ) and lGupta and Vargal (|l993f ). 



Important analogues of the beta funct ion integral are given in the following theorems. 
Theorem 13.91 is discussed bv lKhatril (Il966h in the real case. 
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Theorem 3.9. //Re G'^i'^, then 

f |x|a-(m-l)/3/2-l|j ^ Xr^'^+^'^Cf (RX-l)(dX) 

/xeq3f„ 

_ ri[a,-K]ri[b,n] 
ri[a + b] 

for Re(a) > (to - l)/3/2 + fci and Re(&) > (m - l)/3/2 - fc,„. ^nd 



Cf(R), (3.23) 



Cf(R), (3.24) 



/" |x|'^-(™-i)/3/2-i|i + xr('^+'')cf (RX)(dX) 

_ rg,[a,Ac]rg,[6,-^] 

/or Re(a) > (to - l)/3/2 - fc™ and Re{b) > (to - l)/3/2 + fci. 
Proof. By Corollarv l3.6i we have for any Z e # 

/" etr{-XZ}|X|'^-(™-i)''/2~iCf (X'^R) |Z|"(dX) 

= r^[a,-«;]Cf(RZ). (3.25) 

Multiplying both sides of ((3?25| by etr{-Z}|Z|^"(™"i)^/2"i and integrating with respect to 
Z we have 

f (f etr{-(I + X)Z}|Z|"+''-(™-i)''/2-i(^2)) 

x|X|a-(™~l)/3/2-l(^^ (X-lR) (dX) 

= Ti[a,~K] f etr{-Z}|Z|''-(™-i)''/2-i^^(p^2)(dZ). (3.26) 

The desired result in p.23p is obtained by using Corollary 13.61 and p.l7p in the left and 
right sides of (|3.26p . respectively. The result in (|3.24p is obtained similarly. □ 

Corollary 3.10. //Re 6^'^, then 

I |x|a-(m-l)/3/2-l|j _ X|''-('"-l)'^/2-lCf (RX-l)(dX) 

/o<x<i 

1 m[a + 0, -kJ 

/or Re(a) > (to - l)/3/2 + /ci and Re(&) > (m - l)/3/2. 

Proof. It is obtained in a similar way to that given for p.l3p . see I Gross and Richards! (Il987l) . 

□ 

Now, taking 6 = (m - l)/3/2 + 1 > (to - l)/3/2, from ([XT^ and CoroUarv IXTUl we have 
the following result. 
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Corollary 3.11. //Re ef^"^, then 

[ |xp-('"-i)^/2-iCf(RX-i)(dX) 

"'o<X<I 

ri[a, -«]rg,[(m - l)P/2 + 1] ^0 
Ti[a + {m-l)(3/2 + l,-K] 

for Re(a) > (to - l)/3/2 + ^nd 

/" |X|'^-(™-i)'3/2-icf(XR)(dX) 
io<x<i 

.i-p/3 r^™ 1^/^/o 1 11 

Cf(R), (3.29) 



/o<x<i 

rl[a,K]r^[(TO-l)/3/2+l]^ 



[a+ (m- l)/?/2 + 1,K 

for Re(a) > (to - l)/3/2 - k„,. 

Similarly, taking a — {m — l)(3/2 + 1 > (to — l)/3/2 — fc™ in p.l3p . we have the following 
result. 



Corollary 3.12. //Re 6^*^, then 

I |I-X|''-("-i)^/2-i(;.^(xR)(dX) 
/o<x<i 



= ^-[(-"^)^/^ + ^'-]^-[^]c.^(R), (3.30) 
rf„[(TO-l)/3/2 + l + 5,K] 

/orRe(6) > (to - l)/3/2. 

4 Hypergeometric functions 

In this section, we study diverse integral properties of hypergeometric functions for normed 
division algebras. First, let us consider the following definition. 

Fix complex numbers ai, . . . , Op and 6i, . . . , 6q, and for all 1 < i < g and 1 < j < m do 
not allow —hi + [j — l)(i/2 to be a nonnegative integer. Then the hypergeometric function 
with one matrix argument pFj^ is defined to be the real-analytic function on 6^ given by 
the series 

,F^^ia„ . . . , a,; 6„ . . . , 5, X) = g ^ ["^j; ' ' ' j^^jl (4.1) 

fc=0 K [Oi\k ■ ■ ■ [bq\K 



Some known properties are, see lGross and Richards! (|l987l Section 6, pp. 803-810): 
Convergence hypergeometric functions. 

1. If p < q then the hypergeometric series (14. ip converges absolutely for all X e 6^. 

2. If p = 9 + 1 then the series (|4.ip converges absolutely for ||X|| = max{|Ai| : i = 
1, . . . , to} < 1, and diverges for | |X| | > 1, where Ai, . . . Am are the i-th eigenvalues of 

Xe6^. 



3. If p > q then the series (14. ip diverges unless it terminates. 

For all X e S^; indeed, for all X e ©m*^- This is characteristic of the general situation 
when p < q. 

o^o'^(X) = E E ^ - E ^ = (4-2) 

fe=0 K k=0 
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If Re(a) > (to - l)/3/2, and | |X| | < 1, 
... 1 



iFo^(a;X) 



r„i[a] JYeqj 

II -xr'^ 



etr{-(I-X)Y}|Y 



a-(m-l)/3/2-l 



(dY) 



(4.3) 



gives the full analytic continuation of iF^(a; •) to any simply-connected domain in S^'^. 
The right side is determined by the principal branch of the argument. The fact that the 
hypergeometric series iF^ has {X G 6^^ : ||X|| < 1} as its domain of convergence is 
characteristic of p+iF^ for all p > 0. 

Let Re(c) > Re(a) + (m - l)/3/2 > (to - and | |X| | < 1. Then 



P+i-Ff+i(ai, 



,ap,a;bi, . . . ,bq,c;X) 
1 



Bm[a, c - a] 7o<Y<i 



F'^(ai---ap;6i---6,;XY) 



q 



X |Y|°~^'"~"'"'"'^^~"'"|I — Y|'^"°~'™~"'"-"'/^^"'"((iY) 
for p = q + 1. In particular, for p ^ 2, we have the Eulcr formula 

2Ff (ai,a;c; X) 



(4.4) 



Sm[a, c - a] 



0<Y<I„ 



iFo^(ai;XY)|Yp-("-i)'5/2-i 

X |I- Y|"~"-(™-i)''/2-i(^Y). 



(4.5) 



for arbitrary ai, Re(c) > Rc(a) + (to - l)/3/2 > (m - l)/3 and ||X|| < 1. 

Remark 4.1. Observe that (|4.4p (and of course (|4.5|) . too) is a consequence of (I3.13p . And 
then the condition over a and c must be Re(c) > Re(a) + (m — l)/3/2 > (to — l)/3 — km- 

Laplace transform of hypergeometric functions. Assume p < Re(a) > (m — l)/3/2 and 
U G 6^. Then 



/ etr{-XZ}pi^f (ai • • • a^; 6i • • • 6,; XV)\Xr^"'-'^f'/^-\dX) 
= IZj-T^M p+iFg''(ai---ap,a;&i---6g;UZ-i). 



(4.6) 



When p < q, the integral in (|4.6p converges absolutely for all Z G When p = q, the 
integral converges absolutely for all Z G e^*^, such that ||(Re(Z))-i|| < 1. 

Remark 4.2. iHerzj ( 19551 p. 485) wro te... "The convergenc e of th e integral (|4.6p requires at 
least Re(a) > (to — l)/3/2... "However, Gross and Richard^ ( 1987 ) based their proof of (|4.6p 
in their Theorem 5.9 which is val id for Re(a) > (to — l)(3/2 — km . The s ame remak should 
be considered in the versions of IConstantQ (|l963l ). iMuirheadI (|l982l ). iRatnaraiah et al 
(l2005bh and lLi and X^ (|200i). 



Similarly, the hypergeometric function of two matrix arguments pFq 
the real-analytic function on given by the series 



(m),,3 



oF, 



(m),/3 



(fll, 



2p;b^,...,b,;X,Y)=J2Yl 



is defined to be 
[ai]g ■ • ■ K]g Cg(X)Cf(Y) 



,=0 . 

Some basic properties of (|4.7p are shown below, see iGross and Richards! (|l989l ). 

Convergence hypergeometric functions with two matrix arguments. 



(4.7) 



15 



1. If p < q then the hypergeometric series (|4.7p converges absolutely for all X and Y S 

6/3 

2. If j> = g + 1 then the series (|4.7p converges absolutely for ||X|| • ||Y|| < 1, and diverges 
for||X||.||Y||>l. 



Also 



Heilf (m) 



pFj^iai, . . . , ap; 6i, . . . , 6,; XHYH*)(dH) 



F!;"'^'Piai,...,ap;b,,...,b,;X,Y) 



q 



In particular 



oFo^(XHYH*)(dH) = /" etr{XHYH*}(dH) 

^ o^^o^"^'^(X,Y). 



(4.8) 
(4.9) 

(4.10) 
(4.11) 



Laplace transform of hypergeometric functions with two matrix arguments. 
Assume p < q, Re(a) > (m - 1)13/2 and U e S{^. Then 



y" ^ etr{-XZ}p^^^('")''^(ai---ap;6i---6,;XU,Y)|Xr-(™-i)''/2-i(^X) 
- IZr^r^Ja] p+iF^(™)''3(ai---ap,a;6i---&,;UZ-\Y). 



(4.12) 

When p < q, the integral in (|4.12p converges absolutely for all Z e $ and Y e When 
p = q, the integral converges absolutely for all Z and Y G such that ||(Re(Z))~^|| • 

|Y|| < 1. Remark similar to\ <4. 21 should be considered for (|4.12p . 

We now propose further integral properties of hypergeometric functions for normed di- 
vision algebras. The first result i s the i nverse Lap lace transformati on. In the real case, 
this re sult was obtained bv iHer j ( 19551 ). IConstantinc ( 19631 . Jaincs, (1964) and Muirhead 
( 1982 . p. 261), and in the complex case by Mathai (Il997l p. 370). L et us firs t consider 
the fo llowing extension of similar results discussed in IConstantind (|l963| ) , see also iMuirhead 
(|l98l p. 253). 

Lemma 4.3. Assume that Z and X G &m'^, U £ and Re(a) > oq. Then 

2m(m-l)/3/2 



(-27ri)™("i-l)/3/2+r 



■ f etr{XZ}|Zr'^Cf(UZ-i)(dZ) 

JZ-ZqG* 1 

.|X|a-(m-l)/3/2-1^0(XU), 



(4.13) 



where Zq G 

Theorem 4.4. Assume that Z and X e 6^*^, U £ and Re{b) > bo. Then 

^2li)Lim-,w.+^ l^^^^ etr{XZ}|Z|-%F,^(ai, . . . , a,;b,, . . . , 6,; UZ-i)(dZ) 

= |X|''-(™-i)'3/2-ipFf^,(ai, . . . , ap; 6i, . . . , 6„ 5; XU), (4.14) 

and ifY e 6f„ 



(27rj)™(™-i)^/2+™ 



etr{XZ}|Z|-'' 

Z-Zo6* 

X pi^,('")''^(ai, . . . , ap; &i, . . . , 6,; UZ-\ Y)(dZ) 
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= |X|''-(™-i)'3/2-ipF,(™)'''(ai, . . . , ap; 61, . . . , b„ b; XU, Y), (4.15) 

where Zq G 'Pfn- 

Proof. Proof of both (I4.14p and (|4.15p follows by expanding the pF/^ and pFq™'''^ functions 
in the integrands and integrating term by term using ()4.13|) . □ 

Theorem 4.5. The iF^ function has the integral representation 



iFf (a,c;X) = --i / etr{XY}|Yr-("-i) 

Dm\a, c - a Jo<Y<i 



?^(nr-^\ = / et^-r^V\IVI"-('"-l)/'/2-l 

l] ^0<Y<I 

|j_ Y|'=-'^-(™-i)'3/2-i(dY), (4.16) 



X 



valid for Re(c) > Re(a) + (m - l)/3/2 > (m - l)/3 - fc„ and all X e 6^-*^. 

Proof. The desired result is obtained by expanding etr{XY} using p.Sp and integrating 
term by term using (|3.13p . □ 

The generalised Kummer and Euler relations are given in the following result. 

Theorem 4.6. 

iFf(a,c;X) = etr{X}ii^f (c - a, c; -X), (4.17) 

/orX e 6^*^. And 

2F^{a,b;c;X) = \I - X^F^ {c - a,b; c; -X{I - X)-') (4.18) 
= |I-X|"-"-''2Ff(c-a,c-6;c;X) (4.19) 



/or||X||<l. 



Remark 4.7. Observe that, for iJ^^(a;X) the condition Re(a) > (m — l)/3/2 over a is 
determined by its integral represe ntation (14.31) . However iF^(a;X) is easily seen to be 
analytic for all a and ||X|| < 1, see Herd (ji955, p. 486). Similarly, the conditions Re(c) > 



Re(a) + (m — 1)/3/2 > (to — l)/3 — fc™ over a and c given in Theorem l4.51 valid for Theorem 
(|4.17p too, are determined by the existence of the integral (|4.16p . However, these conditions 
can be extended to other possible values if we use the inverse Laplace transformation to 
define iF^ {a, c;X). In this cas e iF^ j a, c; X ) is valid for the arbitrary complex a, Re(c) > 
(to - l)/3/2 and X G 6^'^, see lHerzj (|l955l p. 487). Also, the conditions Re(c) > Re(a) 



(to - l)/3/2 > (to - l)/3 - km over a and c for and Theorem SUgH]) and (|iTg)) are 
determined by the absolutely convergence of the integral (|4.5p . Again, these conditions about 
a and c can be extended to other possible values using the inver se Laplace t ransformation and 
the results for iFf (a,c;X) obtained as described before, see iHerd (I1955I p. 489). Finally, 
let us take into account that, for any analysis if the integral representation of iF^ {a,c;X) 
or 2Fi{a,b;c;X) is not used explicitly, then the extended conditions for a and c could be 
considered. 

Theorem 4.8. Let X e -C^,.„ and H = (H1IH2) £ U'^(to), Hi G V^_„. Then 

oFf{(3n/2;p^XX*/4) = [ etr(/3XHi)(dH) (4.20) 



[ etr(/3XHi)(dHi) (4.21) 
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Proof. The proof is analogous to t hat given in the r eal case for iMuirhead (Il98l Theorem 
7.4.1) and in the quaternion case Alternative proofs can be established 

in an analogous form to those given by James ( 196 1[ ) and Herz ( 19551 p. 494 -495). For (|4.2ip 
it might be necessary to consider Lemma 9.5.3, p. 397 in iMuirhead (Il982l) . □ 



On the basis of Theorem 13. 4i we now discuss diverse integral properties of generalised 
hypergeometric functions, which contain as particular cases many of the results established 
above. 

Theorem 4.9. Assume p < q and Re(a) > (m — l)/3/2 — km and U G 6^. Then for 
■d = f{z)z''"'+^-'^dz < oo, 

f /(trXZ)pF^'(ai • • • Op; 6i • • • 6,; XU)|Xr-("-i)'3/2-i(dX) 



Hzrrf>] g^Nf---[«.]g[<5(uz;;).,_ 



[bif. ■ ■ ■ [b,f. n 



k]k\ 



(4.22) 



When p < q, the integral in Iji4-^2\ l converges absolutely for all X G When p = q, the 
integral converges absolutely for all Z G &m'^ , such that ||(Re(Z))^"'^|| < 1. 
Similarly, let p < q, Re(a) > (m — l)/3/2 — km and 

V e&t- Then 



I /(tr XZ)pF|'")''3(ai • ■ • flp; 6i • • • 6,; XU, Y)|X|" 



-(t?x-1)/3/2-1 



(dX) 



= |z|-rf„[a]^^ 



a,t---[a,t[a]P (UZ-i)Cf (Y) 



,=0 « [bit---[bX T[am + k]k\ C^(I) 



(4.23) 



When p < q, the integral in ^.231 1 converges absolutely for all Z E ^ and Y G &m- When 
p — q, the integral converges absolutely for all Z and Y G 6™*^, such that | |(Re(Z))~^ 1 1 • 
I|Y||<1. 

Observe that if /(trXZ) etr{-XZ} in Theorem HJl then we obtain (gH) and (|4T^ . 

Theorem 4.10. Assume p < q and Re(a) > {m — 1)13/2 — km and U G &m- Then for 
7 = JJ^ /(z)z'^'"+'=-idz < oo, 

/ /(trXZ)pFf (ai ■■■apM- • • 6,; X-iU)|Xr-('"-i)'3/2-i(dX) 



= |z|-r^^ja] 



k=o . [fol]^••N^ha+(™-l)/3/2 + l]^ 



x.^^.,. (4.24) 



T[am + k]k 



When p < q, the integral in I14-24\ I converges absolutely for all Ti G When p ^ q, the 
integral converges absolutely for all Z G &m^, such that ||(Rc(Z))^^|| < 1. 
Similarly, let p < q, Re(a) > (m — l)/3/2 — km and U G &m- Then 

f /(tr XZ)pF^(")'^(ai • • • a^; 5i • • • b,; X-^U, Y)\Xr^^-''^^^'-\dX) 



|z|-<^r^[a]5^5] 



(-l)^Mg...Mg[a]g 



/3 



C^'(UZ)C;?(Y) 
X ii-^ ^ '^^ -7. (4.25) 



r[am + fc]A:! Cf(I) 
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When p < q, the integral in i4-^5\ ) converges absolutely for all Z € ^ and Y S 6^. When 
p = q, the integral converges absolutely for all Z and Y £ Gfii^, such that ||(Re(Z))~"'^|| • 
l|Y|| < 1. 

Now, let us define /(trXZ) — etr{— XZ} in Tlieorem l4.101 then we obtain: 

/ etr{~XZ}pFl^{ai ■■■ap;bi--- b,; X-^U)\X\--^^-^^f'/^-\dX) 

= IZr'^r^H pi^f+i(ai---ap;6i---65,-a+(m-l)/3/2 + l;-UZ). (4.26) 

and 

/ etr{-XZ}pF^(™)''3(^^ • • • a^; 6i • • • 6,; X-^U, Y)\Xr^"'-'^f'/^-\dX) 

= |Zr"rf„[a] pF^f/{ai ---ap-M--- -a + (m - 1)13/2 + 1; -UZ, Y), (4.27) 

where for both Re(a) > (m — l)/3/2 + fci. 

Similar results to (|4.6p and (|4.12p or (|4.26p and (I4.27P can be obtained from Corollaries 
1X71 and 

Now, we propose the incomplete gamma and beta functions for normed division algebras. 
Theorem 4.11. Let A G G^f- and r2 e Then 

f etr{-AX}|X|"-(™-i)^/2-i(^X) 
"'o<x<r2 

= Bi[a, (to - l)/3/2 + l]\nWFf{a; a + {m - l)/3/2 + 1; -OA), (4.28) 
for Re(a) > (to - l)/3/2 - km- And, let <S <I, then 

f |Yp-(™-i)/3/2-i|i _ Y|''-(™-i"'/2-i(dY) = Bi[a, (to - l)/3/2 + 1] 

^0 



0<Y<S 

X |H|'^2i^f(a,-fo+(™-l)/3/2 + l;a+(TO-l)/3/2 + l;H), (4.29) 
for Re(a) > (to - l)/3/2 - k„i and Re{b) > (to - l)/?/2. 

Proof. For (|4.28p . let us make the transformation X = rZ^/^Rrj^^^ and by applying Lemma 
HHwe have, (dX) = iril(™+i)'3/2+i(dR), with < R < I. Then, expanding etr{-AX} as 
a series of Jack polynomials and integrating term by term using Corollarv l3.11l the desired 
result is obtained. Similarly, (I4.29P is proved by making the transformation Y = H^/^RH^^^ 
from where, applying the Lemma [2.41 we obtain that (dX) ~ |H|'^™+^^'^/^+^(c?R), with < 
R < I, expanding |I - XH|^-("-i'^/2-i = iFl^i-b + (m - l)/3/2 + 1; XH) and integrating 
term by term using Corollarv l3.11l □ 

Theorem 4.12. Let A e &m'^ and $7 6$. If r = a ~ (to — l)/3/2 ~ 1 is a positive integer, 
then 



f etr{-AX}|Xp-("-i)''/2-i(^x) 
Jx>o 



= r^Ja]|A|-" etr{-Af7} ^ ^ (4.30) 

/or Re(a) > (m — l)/3/2 + fci and ^* denotes summation over those partitions k = 
(fci, . . . , km) of k with fci < r. 
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Proof. Consider the transformation X = Jl^/^(I+R)f2-'^/^ and applying Lemma we have, 
(dX) = |n|('"+i)'5/2+i(dR), with R > 0. Noting that |I + R| = |R||I + R-i| and expanding 
|I + R-i|a-("i-i)/3/2-i terms of Jack polynomials, assuming that r = a — {m — l)/3/2 — 1 
is a positive integer we obtain 

\l + ^-i\a-(m-i)i3/2-i ^ ii^(f(-a + (m- l)/3/2 + l;-R-i) 

[-a+ {rn - l)/3/2 + (riA) 

k=0 K 

because [—a + (m — l)/3/2 + = is any part of k that is greater than r. The desired 
result is obtained by integrating term by term using Corollarv l3.6l □ 

We end this section with a some general results, which are useful in a variety of situations, 
which enable us to transform the density function of a matrix X G to the density function 
of its eigenvalues. 

Theorem 4.13. Let X G 'P^j be a random matrix with density function /(X). Then the 
joint density function of the eigenvalues Ai, . . . , Am of X is 

■n(^.-A,)M /(HLH*)(dH) (4.31) 



i<j 



where L = diag(Ai, . . . , Am), Ai > • • ■ > Am > 0, g is defined in Lemma \KR and (dH) is the 
normalised Haar measure. 



Proof. The proof follows immediately from Lemma 12.61 □ 

5 Invariant polynomials 

In this section, we extend many of the properties of a class of homogeneous polynomials for 
normed division algebras of degrees k and t in the elements of matr ices X and Y G 
respectively, see Davis (1979), Davis (1980, ). IChikuse (1980) and Chikuse and Davis! (|l98 



these are denoted as C^''1'''^(X,Y). These homo geneous polynomials are invariant under 
the simultaneous transformations 

X ^ U*XU, Y ^ U*YU, H G U'^im). 

The most important relationship of these polynomials is 

Cf (AH*XH)Cf (BH*YH)(dH) 



L 



HGllf (n 

_ cf'---(A,B)cif'-'-(X,Y) 

where (dH) is the normalised Haar measure and Cf , and are Jack polynomials 
indexed by ordered partitions k, t and 4> of nonnegative integers fc, t and / = k + t, 
respectively, into not more than m parts, (f) G k.t denotes the irreducible representation of 
GL{m,^) indexed by 2(/) that occurs in the decomposition o f the Kroneck er pr oduct 2k ^2 t 
of the irreducible representations indexed by 2k and 2t, see iDavis (Il979ll and lDavid (Il980l) . 



In a similar way to the case of Jack polynomials, let A = A* A and B = B*B. We opt for 
convenience of notation rather than strict adherence to rigor, and write cF'''''^(XA, YB) 
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or CF1''''^(AX,BY) rather than cj?'''''^(AXA*, BYB*), even though XA, YB, AX, or 
BY need not he in 

Some of the elementary properties and results on invariant polynomials are extended 
below: 

Elementary properties of ' . 

Let X and Y e 6^^, then 

cf-^(X,X)=<l''-C,^(X), where = ■ (5-2) 



Cf^'^(X,Y) = 



-Cf(X), forY = I; 



d> !cf(Y), forX^I. 



(5.3) 



cJfl"'°(X,Y) =Cf(X), and cJfl°'"(X, Y) = (Y). (5.4) 
Cf (X)C,^(Y) ^ Yl ^'''"Cf ''•^(X, Y), (5.5) 



therefore, 



n[/3]K,r^[/3]K,- 



(trX)'=(trY)*= 5] 0f'''^cJfl'^'^(X,Y). (5.6) 
From dO]) and (l5?5]) 

Cf(X)Cf(X)= ^ (0P'^'^)'c^(X). (5.7) 

For constant a and 6 

cJfl"'"(aX,6X) = cjf' "^"(X, Y). (5.8) 

The next expansion can be used to derive several useful results of invariant polynomials. 
From ([5?T|) . (fSTS]) and (|4?2l) we obtain 



/ etr{AH*XH + BH*YH}(fiH) 



where 



^ Cf-'-(A,B)Cf-'-(X,Y) 

,.u,^/3.TN ' ^^-^^ 



00 00 CX3 



E=EEEEE- 

H,T;4> k = t=0 K T 0gK.T 



From (|5.9|1 we obtain, see iDfaz-Garcial (|2008l ). 

f c'fl'"'^ (AH*XH,BH*YH)(dH) 

JHeU'^(m) 

c[fl'='^(A,B)cf«^^(X,Y) 



(5.10) 
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in particular using ([5 



/• Ml "■7A*A,B)C^(X) 
/ cfi"'"(A*H*XHA,B)(dH) = ^ ^ — ' i ^\ i 



analogously 

/ cf''='^(A,B*H*YHB)(dH) = ^^ ' (5.12) 

Laplace transform. 

For all A and B G ef„, Z e * 

/ etr{-XZ}|X|°-('"-i)'3/2-i(^[/3lK.r^^^^ BX)(dX) 

= rfja>]|Zr'^cf'''^(AZ-\BZ-i). (5.13) 
valid for Re(a) > (m — l)/3/2 + (fc + In particular 

/ etr{-XZ}|X|"-(™-i)^'/2-iC^fl'''^(AXA*,B)(dX) 

= r^[a,«:]|Z|-''cf''-^(AZ-iA*,B), (5.14) 

where Re(a) > (m - l)/3/2 + /ci. And 

/ etr{-XZ}|Xp-("-i)''/2-i(^[/3]«^,r^A, BXB*)(dX) 

= r^[a,r]|Zr'^cf'''^(A,BZ-iB*), (5.15) 

with Re(a) > (to - l)/3/2 + h. 

Similarly, for aU A and B G e^, Z e *, 

/ etr{-XZ}|Xp-('"-i)'3/2-icJ'l''^''(AX-\ BX-i)(dX) 

= r^[a,-0]|Z|-"Cfl'''^(AZ,BZ), (5.16) 
where Re(a) > (m — l)/9/2 — (fc + t)m. In particular 

/ etr{-XZ}|X|°-(™-i'^/2-icJf''''7AX-iA*,B)(dX) 

= Ti[a,-K]\Zrcf-'^iAZA*,B), (5.17) 

with Re(a) > (to - 1)^/2 - fc^. And 

/ etr{-XZ}|X|»-(™-i)'3/2-icf l^'^A, BX-iB*)(dX) 

= r^[a,-T]|Zr'^cJfl«^^(A,BZB*), (5.18) 

valid for Re(a) > (m - l)/3/2 - i^. 
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Inverse Laplace transform. 

Assume that Z and X G 6^;*^, A and B e 6^, and Re{b) > bo. Then 



(27rz)™("-i)/3/2+ 



f etr{XZ}|Zr''Cjfl'''''(AZ-\BZ-i)(dZ) 
Jz-Zoe* 

^ |X|&-("-i)'3/2-ic|fi"'"(AX,BX), (5.19) 



and 



(•27ri)'"(™-i)/3/2+ 



/" etr{XZ}|Z|-''Cjfl'''^(AZ,BZ)(dZ) 

^ |X|^-(™-i)'3/2-ic|fl'''"(AX-\BX-i). (5.20) 



Similar expressions are obtained for c|f' '''''(AZ-\ B) and cjf^'"^^(A, BZ"^) from (OQ]) : 
and for cJf]"^"(AZ, B) and cJf]"^"(A, BZ) from (IQUll . 

Beta type I integrals. 

For all A and B e Sf^*^ and Re(6) > (m - l)/3/2, 

/ |Xp-("-i)/5/2-i|i_X|''-("-i)^/2-icP"'"(AX,BX)(dX) 
-'o<x<i 

= W^^r^(A,B), (5.21) 

valid for Re(a) > (m — l)/3/2 — (fc + i)m. In particular 

/ |X|--(™-iW2-i|i_x|''-(™-iW2-ic[fl'''"(AXA*,B)(dX) 
Jo<x<i 



/0<X<I 

_ r^[a,«;]r^[&] 



Tmia + b, k] 

with Re(a) > (m - l)/3/2 - fc,„. And 



C7'f'"'"(AA*,B), (5.22) 



/ |X|«-(™-i)/3/2-i|i_x|^-("^i)'5/2-i^[/3]-,r(A^gXB*)(dX) 
Jo<x<i 

^ r|klHtMc'fl«^"(A,BB*). (5.23) 

where Re(a) > (m — l)/3/2 — 1,„. Another particular integral given in the real case bv I David 
(Il979l ) is 

/ |X|-(™-iW2-i|i _ x|f-("-i)/3/2-ic[fl---(x,I _ X*)(dX) 
Jo<x<i 

= %r%^<'"^^^.'(i)' (5-24) 

ri;[a + 6, 0J 

valid for Re(a) > (m - l)/3/2 - fc„j and Re(6) > (m - l)/3/2 - fc^- 
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Analogously, for all A and B G Ql^f and Rc(6) > (to - l)/3/2, 

/ |X|"-('"-iW2-i|i_x|''-(™-i)'3/2-icJfl"'"(AX-\BX-i)(rfX) 
Jo<x<i 

Tm[a + 6, - 

where Re(a) > (m — l)/3/2 + (fc + In particular 

^0<X<I 



^^^^ -C'r"^{A,^), (5.25) 



/0<X<I 

re,[a,-A«]rf„[6] 



[a + 6, -k] 

valid for Re(a) > (m - l)/3/2 + fci. And 



C7''''"^"(AA*,B), (5.26) 



I |x|«-(™-i)/^/2-i|i _ x|''-(™-i)''/2-iC7jfl'='"(A, BX-iB*)(dX) 
^o<x<i 

d'^l"'"(A,BB*), (5.27) 



/o<x<i 

rf„[a + 6,-r] 

with Re(a) > (m - l)/3/2 + ii. 

Now, taking b= {m- l)/3/2 + 1 > (m - l)/3/2 in (|OT|l and we have the following 

results. 

For all A and B e e^*^, 

f |xr-(™-i)'3/2-icJfl'^'"(AX,BX)(dX) 

"'o<X<I 

= rg.K^]rg.[(--lW2 + l]w.K.(^^B), (5.28) 
r^[a+(m-l)/3/2 + l,0] ^ ^ ^ 

valid for Re(a) > (m - l)/3/2 - (fc + *)„. And, 



Jo<x<i 



r^L[a,-'^]r^^[(--lW2 + ll^[,,..(^^B), (5.29) 



ri[a+{m-l)P/2 + l,- 
where Re(a) > (m - l)/3/2 + (fc + t)i. 

Beta type II integrals. 

For aU A and B e 6^^, 

f |x|«-(™-i)'5/2-i|i + xr c[fl"'"(AX, BX) 

= liMHil^cf 1-^(A, B), (5.30) 
T?ri [a + b] 

with Re(a) > (to - l)/3/2 - (fc + t)„ and Re(6) > (to - l)/3/2 + (fc + t)i. In particular 

f |x|a-(m-i)/3/2-i|i ^ xr('^+^)c|fl'^'"(AXA*,B) 
Jxe?pl 

^ rg,[a, A.]rgJb, -^] g[/3]^,^(AA*,B), (5.31) 
T'mia + b] 
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such that Re(a) > (m - l)/3/2 - km and Re(6) > (m - l)/3/2 + fci. And 
/ |X|«-(™-i)'5/2-i|I + X|-("+'')c[fl"'"(A,BXB*) 

i(7[fl'''^(A,BB*), (5.32) 



ri[a,r]Ti[b,-r] 



Ti[a + b] 

vahd for Re(a) > (m - l)/3/2 - t„ and Re(&) > (m - l)/?/2 + ti. 
In a similar way, for all A and B G S^^, 

f |xr-(—i)'5/2-i|I + X|-(«+*')cJ^l'='"(AX-\BX-i) 

^ ri[a,-ct>]Ti[b,^] ^[^]^ 
Ti[a + b] 



^^^^ Cr'^'^(A,B) (5.33) 



where Re(a) > (m - l)/3/2 + {k + t)i and Re(6) > (m - l)/3/2 - (A; + 1)^. In particular 

/ |Xr-('"-i)'3/2-i|i^x|-('^+'')c|fi"'"(AX-iA*,B) 

with Re(a) > (m - l)/3/2 + fci and Re(6) > (m - l)/3/2 - And 
f |-j^|a-(m-i)/3/2-in ^ xr(»+^)cf "'"(A, BX-^B*) 



'xeqjg, 

_ r^[a,-T]r^[6,r] 



CfJ'''^(A,BB*), (5.35) 



such that Re(a) > (m - l)/3/2 + ti and Re(6) > (m - l)/?/2 - tm- 

Incomplete gamma and beta functions. 

First consider the following results 
For all A and B e and < H < I, 

r |X|«-(™-i)'5/2-icf "'"(AX, BX)(dX) 

= ^">'^]^-[(--^)^/^ + ^]|»rcr--(AB,BB). (5.36) 

r^[a+(m-l)/3/2+l,<^] ' I ^ ^ ' ^ 

valid for Re(a) > (m - l)/3/2 - (/c + t)m- And 

/■ iXI^-^^-i^'^/^-ijoJJl^.r g^^^X^ 

^0<X<S 

^ r(^K^]rg,[(rn-i)/?/2+i] i,,.^^ 

r^[a+(m-l)/3/2+l,K] ' ' ^ ^ 
valid for Re(a) > (to - l)/3/2 - km- 
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The next result is obtained immediately, expanding etr{— XA} in terms of Jack polyno- 
mials, making use of the property (|5.5p and integrating term by term using (|5.36[) . Thus, 
for all A and B e ef„ and G *, 



J etr{-XA}|X|"-('"-i)'^/2-iC^(BX)(dX) 

r^[a + (m-l)/3/2 + l] t^o.^^^.r k\[a + {m - l)p/2 + 1 



valid for Re(a) > (m - l)/3/2 + (fc + 

Similarly, we expand |I — X|''~(™~^)''/^~^ in terms of Jack polynomials,make use of the 
property (|5.5p and integrate term by term using (|5.36p . For all A S S^;*^ and < H < I, 

j-j^|a-(m-l)^/2-l|j _ -j^|6-(m-l)/3/2-l^/3/^-j^w^-j^x ^ r^[a]r^[(m - l)/3/2 + 1] 

o<x<H rK n ) r'L[a + {m~l)(3/2 + l] 

^ Ij^iay^ M+(^-lW2 + l]X'"^Cy-(»,AO) 

fro.;^... fc![a+(m-l)/3/2 + l]^ 

valid for Re(a) > (m - l)/3/2 - (fc + 



6 Application 

As an application, in this section we found the joint density eigenvalue of the central Wishart 
distribution for normed di vision algebras, and also deri v ed the largest and smallest eigenvalue 
distributions. First, from iDiaz-Garcia and Gutierrez (|2009bh let us consider the following 
definitions. 

Definition 6.1. Let X E Cmn be a random matrix. Then X is said to have a matrix 
variate normal distribution X ~ JV^y.„^{fi, S, 0), of mean fi and Cov(vecX-^) = (g) S, if 
its density function is given by 

^ etr|-^S-i(X-/x)*0-i(X-/x)l . 

(27r/3-i)'^""/^|SP"/2|0|/3W2 i 2 ' y i-Jj 

Also 

Definition 6.2. Let X G £^^„ with distribution X - J\ff^,,^{0, S, I„) and define S = X*X, 
then S is said to have a central Wishart distribution S ^ ("-j ^) with n degrees of freedom 
and parameter S. Moreover, its density function is given by 

|S|/3(n-m+l)/2-l ctr{ " ^ S/2} 



(2/3-1)'^'""/' r^[/3n/2]|S|/3"/2 
with n > (to — 

Therefore, from (|4.3ip and (14. lip the joint density of the eigenvalues, Ai > • • • > Am > 0, 
of S is 



(2/3-1)'"""/' r^[/3n/2]r^[/3TO/2]|S|/3»/2 f^-^ 

[](A,-Aj)Vo(-/3S-V2,L) 



X 

i<j 
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where L = diag(Ai, . . . , A„i). 

In addition, as an immediate consequence of Theorems 14.111 and 14.121 we obtain the 
following result. 

Theorem 6.3. Let S ~ W^(n, S) and n e then 

r^[(m-l)/3/2 + l] |f2|/3"/2 



p(s < n) 



X iFf (/3n/2; (n + m - l)/3/2 + 1; -/3f2S-V2), (6.1) 



valid for Iic(n) > (m — l)/3 — 2fc,„. And if r ~ (n — m + l)/3/2 — 1 is a positive integer, then 
Pin > S) = etr{-/3»S- V2} f] ^ * ^'(^"^f ' ^2) 3) 

fe=0 K 

/orRe(n) > (m— l)/3+2fci a7^rfX]K denotes summation over those partitions k — [ki, . . . ,km) 
of k with ki < r. 

Observing that if Amax and Amin are the largest and smallest eigenvalues of S, respectively, 
then the inequalities Amax < x and > y are equivalent to S < and S > yl, 

respectively and the following result is obtained. 

Corollary 6.4. Assume that S ^ Wm(?^, S) and x > 0. Then 

rilim- 1)13/2+1] a;/3mn/2 



-P(A,nax < x) 



(2/3-1)"""/' Ti[{n + m- 1)13/2 + 1] ISI^'"/^ 

X iFf{Pn/2; (n + to - l)(3/2 + 1; -(3x1:-^ /2), (6.3) 



valid for Re(n) > (to — l)f3 — 2km- And if r — (71 — to + l)(3/2 ~ 1 is a positive integer and 
y > 0, then 

P(Ami„ <y) = l' etr{-/3j;S- 72} E E (6-4) 

fc=0 K 

/or Re(n) > (m— l)/3+2fci and ^* denotes summation over those partitions n ~ (fci, . . . , km) 
of k with ki < r. 

As a numerical example we plot the distribution function of Xmax on Figure 1 and the 
distribution function of Xmin on Figure 2. First note that applying the generalised Kummer 
relation (|4.17p in (|6.3p we obtain 

^ rgJ(TO - l)/3/2 + 1] etr{-/3xS-V2} o^""'"/^ 

lAmax < X) ^^^^^^,mn/2 ^ ^ _ ^-^^^^ ^ |S|/3n/2 

X li^f ((to - l)/3/2 + 1; (n + TO - l)/3/2 + 1; f3x'E-'^/2). 
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The smallest eigenvalue distribution of real central WIsharl distribution The smallest eigenvalue distribution ot complex central WIshart distribution 




2 4 6 8 10 12 2 4 6 8 10 12 



Figure 2: Distribution functions of Amin of (7, diag(l, 2)), /3 = 1,2,4 and 8. 
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